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i Recent results on identified hadrons from the PHENIX experiment in Au+Au coUi- 
CN I sions at mid-rapidity at -^/snn = 200 GeV are presented. The centrality dependence of 
transverse momentum distributions and particle ratios for identified charged hadrons are 
■ studied. The transverse fiow velocity and freeze-out temperature are extracted from pt 
spectra within the framework of a hydrodynamic collective fiow model. Two-particle HBT 
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O '. correlations for charged pions are measured in different centrality selections for a broad 
^ I range of transverse momentum of the pair. Results on elliptic flow measurements with 
^ ■ respect to the reaction plane for identifled particles are also presented. 

o ■ 

53 ; 1. INTRODUCTION 

CJ I The physics motivation of the ultra-relativistic heavy-ion program at the Relativistic 
^ ■ Heavy Ion Collider (RHIC) is to study nuclear matter at extremely high temperature 
and energy density with the hope to reach a new form of matter called the quark gluon 
plasma (QGP). Among the various probes of the QGP state, hadrons carry important 
rS I information about the collision dynamics along with the spatial and temporal evolution 
I of the system from the early stage of the collisions to the flnal state interactions. 

For studies of hadron physics at RHIC, the PHENIX experiment [ |I[ demonstrates 
good capability for particle identiflcation (PID) for both charged hadrons (vr^, K^, p, 
p, d and d) and neutral pions over a broad momentum range. The charged hadrons 
can be identifled with time-of-flight measurements in two different detectors: (1) a high 
resolution Time-of-Flight wall (TOP) and, (2) an electro-magnetic calorimeter (EMC), in 
conjunction with the tracking system in the PHENIX central arm spectrometers and the 
beam counter, which provides the start timing and the event vertex determination. The 
tracking system in the central arm consists of drift chambers (DC), three layers of pad 
chambers (PC), and time expansion chambers (TEC). The PHENIX central arms cover 
|?7| < 0.35 in pseudo-rapidity, and cover 7r/4 with the TOP and 371/4 by EMC in azimuth. 
The tt/K and K/p separation can be achieved up to 2 and 4 GeV/c in pt-, respectively, 
using the TOP detector, which has a 120 ps timing resolution. The EMC offers a larger 
azimuthal coverage for PID, but suffers from a timing resolution of only 500 ps. Neutral 
pions are identifled with the EMC via the 7r° — ^ 77 decay channel up to 10 GeV/c in px 
using the full statistics taken during Run II at RHIC in Au-|-Au collisions [ 0] . 



*for the full PHENIX Collaboration author list and acknowledgements, see Appendix "Collaborations' 
of this volume. 
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Figure 1. Transverse momentum distributions for pions (circles), kaons (squares) and p, p 
(triangles) in the 0-5% most central events (upper panels) and 60-91.4% most peripheral 
events (lower panels) at y/s^^ = 200 GeV in Au+Au collisions. The left panels show pos- 
itive particles and the right panels show negative particles . The error bars are statistical 
only. 



2. IDENTIFIED SINGLE PARTICLE SPECTRA 

We have measured the transverse momentum distributions for vr^, K^, p and p at mid- 
rapidity in Au+Au collisions at ^/snn = 200 GeV over a broad momentum range over 
various centrality selections. To identify the charged particles, the high resolution TOF 
counter is used in this analysis. We use about 4 million minimum bias events. The data 
are classified into 11 centrality bins expressed in percent of the total inelastic cross section. 
The spectra for each particle species are corrected for geometrical acceptance, decay in 
flight, multiple scattering, and tracking efficiency using a single particle Monte Carlo 
simulation. A multiplicity-dependent track reconstruction efficiency is also determined 
and applied by embedding simulated tracks into real events. 

The upper two panels in Figure ^ show the px distributions for identified hadrons in the 
most central events (0-5%) in 200 GeV Au+Au collisions for positive (left) and negative 
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Figure 2. Mean transverse momentum for identified charged hadrons as a function of 
Npart for protons and anti-protons (triangles), kaons (squares) and pions (circles). The 
left panel shows positive particles and right panel shows negative particles. The open 
symbols indicate the data for 130 GeV [ ^ while the filled symbols indicate the data for 
200 GeV in Au+Au collisions. 



(right) particles. The lower two panels show the most peripheral events (60-91.4%). In 
each panel, the data are presented up to 1.8 GeV/c for charged pions and kaons, and 3.8 
GeV/c for p and p. In the low pt region of the most central events, the data indicate 
that the inverse slope increases with the particle mass. Also, the shape of the p and p 
spectra have a shoulder-arm shape while the pion spectra have a concave shape. On the 
other hand, in the most peripheral events, the spectra are almost parallel to each other. 
This mass dependence of the slopes and shapes of the spectra for protons in the central 
events can be explained by a radial flow picture. At around 2.0 GeV/c in central events, 
the proton yield becomes comparable to the pion yield. A similar behavior is also seen 
for negatively charged particles. 

In order to quantify the observed mass dependence of the slopes, the mean transverse 
momenta (pt) are extracted for 11 centrality selections and for each particle species. The 
Pt spectra are extrapolated to below and above the measured range using a power law 
function for pions, an rrij- exponential function for kaons, and a Boltzmann function for 
p and p. In Figure H, the centrality dependence of {px) for identified charged hadrons are 
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Figure 3. Freeze-out tempera- 
ture (top) and transverse flow 
velocity (bottom), func- 
tion of Npart, extracted from 
the hydrodynamical collective 
flow model fit to the single par- 
ticle spectra at ^/sNN = 200 
GeV in Au+Au. 



shown together with the 130 GeV data points [ In both the 200 GeV and 130 GeV 
data, (pt) for all particle species increases from the most peripheral to the most central 
events and also increases with particle mass. The dependence of (pt) on particle mass 
suggests the existence of a collective hydrodynamical expansion. The dependence on the 
number of participant nucleons (Npart) calculated using a Glauber model [ ^ may be due 
to an increasing radial expansion from peripheral to central events. 

Using these pt distributions, one can characterize the collision system with a small 
set of hydrodynamic quantities based on the expanding source model [ transverse 
flow velocity (Pt) and freeze-out temperature (Tfo). By using this model, which assumes 
boost invariance and a linear velocity profile, a simultaneous fit has been performed for 
the measured spectra in 200 GeV Au+Au in eleven different centrality selections [|^. For 
the 0-5% central events, an expansion flow velocity on the surface, Pt, of 0.7 ± 0.2 (syst.) 
and a common freeze-out temperature, Tfo, of 110 ± 23 (syst.) MeV is extracted. The fit 
results of all particles within each event centrality are shown in Figure ^. The top panel 
is for Tfo and the bottom panel is for Pt, both plotted as a function of Npart- Within 
the systematic uncertainties, the expansion parameters Tfo and Pt decrease and increase, 
respectively, with the number of participants, saturating at mid-centrality. 



3. TWO-PARTICLE HBT CORRELATIONS 

Another way to extract information about the particle emitting source created in rela- 
tivistic heavy ion collisions is with two-particle HBT correlations. The HBT measurement 
is sensitive to the space-time evolution and duration time of the system at the freeze-out 
stage. Using the large data sample of pion pairs collected in Run II, we extend our HBT 
measurements up to 1.2 GeV/c in the transverse momentum of the pairs (fcy) for pi- 
ous [ 1^ . In this analysis the Bertsch-Pratt parameterization is employed in a longitudinal 
co-moving system (LCMS), where the three-dimensional Gaussian radius parameters are 
Rside, Rout and Riong- Thcsc radius parameters are studied as a function of fc^ and cen- 
trality. We use 50 million minimum bias events with charged pions identified using the 
time-of-fiight measurement in the EMC. A full Coulomb correction is applied assuming a 
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Gaussian source and no pairs coming from resonance decays. 

As shown in reference [ ^, we have observed a clear kT dependence for all Bertsch- 
Pratt radius parameters for pions in 200 GeV Au+Au collisions. This kj- dependence 
of the radius parameters can be explained by a space-momentum correlation effect due 
to the expansion of the system [ |^ . For the study of the duration time of the system, 
the ratio Rout/Rside can be used. In Figure ^, the ratio Rout/Rside as a function of kx 
(left) and Npart for pion pairs is shown. The ratio does not change as a function of both 
/ct and centrality within the experimental uncertainties. In contrast to the success of 
the hydrodynamic model for pt distributions in Section 2 of this paper, the HBT radii 
from the hydrodynamic calculation agree only qualitatively with the data with significant 
quantitative discrepancies [ |11|]. At this time, no minimum could be found for the 
simultaneous hydrodynamic fit of the particle spectra and the k^ dependence of HBT 
radius parameters for pions [ O . 
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Figure 4. The ratio Rout/Rside for pions as a function of kj^ (left) and Npart (right) with 
statistical error bars (200 GeV data only) and a systematic error band. In the left figure, 
results from PHENIX and STAR at 130 GeV are also shown [ |, 



4. DEUTERON AND ANTI-DEUTERON SPECTRUM 



Deuteron and anti-deuteron measurements provide complementary information to the 
HBT measurements because the coalescence coefficient B2 [ [12| is considered to be pro- 
portional to the inverse of the volume of the system, and the dynamical evolution of 
the system can be studied with the pt dependence of B2, similar to two-particle HBT 
measurements. 

Figure |^ shows the px distributions of d (squares) and d (triangles) for minimum bias 
events in 200 GeV Au-|-Au collisions, identified with the high resolution TOF detector. 
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Figure 5. px distributions for deuterons 
(squares) and anti-deuterons (triangles). 
The error bars are statistical only. The lines 
are fits using an ttlt exponential function. 



Figure 6. px dependence of the co- 
alescence parameter, B2, for deuterons 
(squares) and anti-deuterons (triangles) at 
mid-rapidity. 



The fitted lines using an ttit exponential function are also shown. We obtain an effective 
temperature of 515 ± 26 (stat.) MeV for deuterons and 488 ± 26 (stat.) MeV for d. 
Figure ^ shows the transverse momentum dependence of B2 for deuterons (square) and 
anti-deuterons (triangle) at mid-rapidity. B2 increases linearly with px up to 3.5 GeV/c. 
The pt dependence of B2 is consistent with the expanding source picture. A similar 
behavior has been observed in the HBT radius parameters for pions [|^. 

5. PARTICLE RATIOS 

Particle ratios provide us with information about the chemical properties of the collision 
system. In Figure the particle ratios of 7i~ /tt^ (left) and / (right) are shown as 
a function oipT for central (0-5%) and peripheral (60-91.4%) Au-|-Au and proton-proton 
collisions at ^JsMN = 200 GeV. A similar plot for the p/p ratio is shown in Figure |^ [ 0]. 
Regardless of the particle species, centrality, or collision system, the particle ratios are 
almost flat as a function of pt within the systematic errors and measured px range, with 
the exception of the p/p ratio for peripheral events, which seems to decrease in the high 
Pt region. The integrated ratios over the measured pt range in the most central events in 
200 GeV Au+Au collisions are: 1.02 ± 0.02 (stat.) ± 0.1 (sys.) for 7r-/7r+, 0.92 ± 0.03 
(stat.) ± 0.1 (sys.) for K-/K+, and 0.70 ± 0.04 (stat.) ± 0.1 (sys.) for p/p. Based on a 



simple statistical thermal model [ |T^|, a baryon chemical potential, /i^, of ~ 30 MeV for 
200 GeV Au+Au central collisions is estimated from the K~/K^ and p/p ratios. 

Particle composition at high pt is also interesting in order to understand baryon pro- 
duction and transport, system evolution, and the interplay between soft processes and 
jet quenching in hard processes. In Figure the p/n and p/ir ratios are shown. For 
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Figure 7. Particle ratios of 7r^/7r+ (left) and K^/K^ (right) as a function of px- The 
data points are for the most central (filled circles) and the most peripheral (open squares) 
in Au+Au and proton-proton (open diamonds) collisions at a/sjvtv = 200 GeV. 



the denominator of the ratio p/n and p/vr, we use two independent measurements with 
different subsystems: (1) px spectra for vr^ up to 2 GeV/c identified using the TOF, 2) 
neutral pions results from 1 GeV/c to 3.8 GeV/c measured using the EMC. Both the 
p/ir and p/n ratios have a clear centrality dependence. The data shows that these ratios 
in central collisions reach unity at a pt of 2~3 GeV/c, while they saturate at around 
Pt of 0.3 - 0.4 in peripheral collisions. This observed behavior in central events may be 
attributed to the composition of two effects including a larger flow effect for protons (and 



p) compared to pions, and a pion suppression effect at high px [ |T^. 



6. ELLIPTIC FLOW FOR IDENTIFIED HADRONS 

In non-central collisions, the initial overlap region of two nuclei is elliptically deformed in 
the transverse plane, resulting in anisotropic pressure gradients. These cause a more rapid 
expansion into the reaction plane than perpendicular to it, resulting in an anisotropy of 
the final px distributions called elliptic flow [ |lT|. Since the event anisotropy is considered 
to be developed at the early stage of collisions, the study of elliptic flow provides a good 
tool to investigate the possible formation of a QGP state. 

Elliptic flow is quantifled by the second harmonic coefficient {V2) of a Fourier expansion 
in the azimuthal distributions of the measured spectrum with respect to the reaction 
plane. We have measured the V2 parameter for identifled particles with respect to the 
reaction plane, which is deflned in the beam counters (|?7| = 3 ~ 4) in 200 GeV Au+Au 
collisions. The details of the analysis method and results are shown in [ 16|. Figure 10 
shows the px dependence of the V2 parameter for identifled particles with respect to the 
reaction plane for minimum bias events. The solid circles are for protons and p, and the 
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Figure 8. p/p ratios as a function of pt for central 0-5% (filled circles), peripheral 60- 
91.4% (open squares) Au+Au and proton-proton (open diamonds) collisions at ^/snn = 
200 GeV. The error bars indicate the statistical errors and the lines indicate the systematic 
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open triangles are for the combined results for pions and kaons. The left panel shows 
negatively charged particles and the right panel shows positively charged particles. The 
solid lines are the results of a hydrodynamic calculation including a first order phase 
transition with a freeze-out temperature of 120 MeV [ 0. The data shows that at lower 
Pt (< 2 GeV/c), the light mesons have a larger V2 parameter compared to protons and 
p. The model calculation agrees very well with the data for all particles in this pt range. 
However, the data seems to deviate from the hydrodynamic calculations above 2 GeV/c 
for mesons and baryons. In addition, the model cannot reproduce the opposite mass 
dependence of the parameter, which is observed above 2 GeV/c. 
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Figure 10. Transverse momentum dependence of V2 with respect to the reaction plane for 
(vr" + K^) and p (left), and for (vr^ + K^) and p (right). The solid lines represent the 
hydrodynamic calculation [ |I7| for pions (upper curve) and protons (lower curve). 



7. SUMMARY 



We present results on identified hadrons in Au+Au collisions at ^JsNN = 200 GeV at 
mid-rapidity over different centrality selections from the PHENIX experiment. The trans- 
verse momentum distributions for vr^, K^, p, p, d, and d are measured and we observe 
a mass dependence of the slopes and shapes of the identified charged spectra in central 
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events. However, they are almost parallel to each other in the most peripheral events. 
Also, the mean transverse momentum for all particle species increases from peripheral to 
central events, and with particle mass. In two-particle HBT correlations for pions, there 
is a clear dependence of all radius parameters. These results on the single particle 
spectra and two-particle HBT correlations are qualitatively consistent with a hydrody- 
namic collective expanding source picture, although there are quantitative discrepancies 
between the model and the HBT results. To characterize the collision system created in 
Au+Au collisions at RHIC energy, the flow velocity and freeze-out temperature have been 
extracted within the framework of a hydrodynamic expanding source model. 

We also present results on the particle ratios 71^ /n^, /K~^, and p/p as a function of 
Pt for different centralities in Au+Au and proton-proton collisions. It is found that the 
ratios are almost constant regardless of particle species, centrality, and collision system. 
The ratios p/ir and p/ir up to 3.8 GeV/c are also measured by using the combined results 
for charged and neutral pions. In central collisions, the ratio reaches unity at a pt of 2~ 
3 GeV/c and saturates at a pr around 0.3 - 0.4 GeV/c in peripheral collisions. 

Finally, we present elliptic flow measurements with respect to the reaction plane for 
identified particles. The data show that at low pr (< 2 GeV/c) the light mesons have 
a larger V2 parameter compared to protons and p. The hydrodynamic model calculation 
agrees very well with the data for all particles up to 2 GeV/c. However, the data seems 
to deviate from the hydrodynamic model above 2 GeV/c for both mesons and baryons. 
The model cannot reproduce the opposite mass dependence of the V2 parameter above 2 
GeV/c. 
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